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The gauge dependence of some fermion bilinear RI/MOM renormalization constants is studied by comparing 
data which have been gauge-fixed in two different realizations of the Landau gauge and in a generic covariant 
gauge. The very good agreement between the various sets of results and the theory indicates that the numerical 
uncertainty induced by the iattice gauge-fixing procedure is below the statistical errors of our data sample which 
is of the order of (1 1.5)%. 



1. Introduction and Strategy 

Non-perturbative renormalization techniques 
(l|^ have become a crucial ingredient in lat- 
tice computations of fundamental QCD pa- 
rameters and hadronic matrix elements. The 
RI/MOM non-perturbative renormalization pro- 
posed in Refs. has been successfully ap- 
plied to compute renormalization constants of 
composite fermion operators in many lattice reg- 
ularizations |^Jl|,^,||j^. The method imposes 
renormalization conditions on conveniently de- 
fined amputated Green functions computed non- 
perturbatively between off-shell quark states at 
large virtuality in a fixed gauge: 

Zr{p\Or\p)\p.=^2^{p\Or\p)o- (1) 

The renormalization constants depend in gen- 
eral on the external states and therefore on the 
gauge in which the renormalization conditions 
have been imposed. For gauge-invariant opera- 
tors the coefficients needed to match in a given 
gauge-invariant scheme has to cancel these de- 
pendences up to higher orders in the continuum 
perturbative expansion and up to discretization 
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errors. 

The existence of both continuum and lattice 
Gribov copies and the numerical noise that they 
can generate (for a recent review see Ref. [|| ) is an 
unsolved problem of the lattice non-perturbative 
gauge fixing. The real concern is, of course, 
the influence that these phenomena may have on 
the values of the renormalization constants, when 
computed using non-gauge- invariant quantities as 
in the RI/MOM scheme. 

Here we discuss the results of a numerical study 
of the systematics induced by the gauge-fixing 
procedure on the RI/MOM determinations of the 
quark field {Z^), the axial- vector {Za) and of the 
scalar density (Zs) renormalization constants (for 
more details see Ref. |^). We have done this 
by imposing the standard lattice Landau gauge 
and another realization of the Landau condition 
(Landaul gauge). Whereas the two realizations 
impose the same gauge-fixing condition in the 
naive continuum limit, they are affected by dif- 
ferent Gribov ambiguities. By comparing the dif- 
ferent sets of results for the renormalization con- 
stants, we have found differences which are negli- 
gible within our statistical errors (of the order of 
(1^1.5)%). 
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An interesting feature of the Landaul gauge is 
that it can be generaUzed to impose a generic 
covariant gauge on the lattice, as proposed in 
Ref. [l0| . By exploiting this opportunity, we have 
performed an exploratory study of the gauge de- 
pendence of off-shell Green functions measured in 
a generic covariant gauge. 

2. Results and Discussion 

We have used a sample of 80 SU(3) gauge 
configurations retrieved from the repository at 
the "Gauge Connection" (http://qcd.nersc.gov/), 
which were generated with the Wilson gluonic ac- 
tion at /3 = 6.0 and y X T = 16^ x 32. By using 
the discretized version (see @ for details) of the 
following functionals: 



H4G]^ J d^xTT[id,A^r] 



H(A,A)[G] 



d^xTv[{d,A<;-Af] 



(2) 



(3) 



(4) 



we have rotated each configuration in the Landau 
gauge, in the Landaul gauge and in the covariant 
gauge. The function A belongs to the Lie alge- 
bra of the group and has been generated with a 
gaussian distribution corresponding to the value 
^ = 8 of the bare gauge fixing parameter. 

Once the quark propagator S{x,0) has been 
computed for each SU(3) configuration and 
gauge and Fourier-transformed, we have deter- 
mined the inverse Euclidean bare-quark propa- 
gator S^^{pa). We have computed the ampu- 
tated Green functions of the local quark bilinears 
Ot{x) = ip{x)Tip{x) in the Fourier space (in what 
follows we adopt the conventions of Ref. pl|) 



Aripa) = S'\pa)Gripa)S-\pa) , 



(5) 



where 



Grip) 



d'^xd^y e 



-ip{x-y) 



{i:{x)Orm{y)) (6) 



and r is a generic Dirac matrix corresponding 
to F = {A, S} = {7^75, 1}. The renormalization 




Figure 1. The Za and Zs renormalization constants 
are shown as a function of (a/i)^ for the Landau 
(crosses) and Landaul (triangles) gauges. The data 
have been slightly displaced in the ar-direction and the 
error bars for one set only have been shown. 



constant Z^^{fia, go), which defines the renormal- 
ized operator O^^ = Z^^Or, is fixed by imposing 
in the chiral limit the renormalization condition 



Z^\^ia)ZT\^m)^T PrAr(pa)|p2^^2 = 1, 



(7) 



where Pp is a suitable projector on the tree- 
level amputated Green function and where the 
RI wave-function renormalization constant Z^ is 
defined as 



Z^l^a) — lim 



dS~^{pa) 



■ (8) 

Even though Oa and Os are gauge- 
independent operators, their matrix elements 
between quark states are gauge dependent. In 
Fig. |l|, the numerical results for Z^, Za and 
Zs calculated in the Landau and in the Landaul 
gauge as a function of the square lattice momenta, 
are shown. Our data are in very good agreement 
with the results reported in the literature [Tl|| . 

The results for all correlators, corresponding 
to Landau and Landaul gauges, coincide within 
the statistical errors. As a consequence we can 
conclude that the data do not show any resid- 
ual effect due to the presence of lattice Gribov's 
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Figure 2. The behaviour of Za and Zs is shown as a function of [a/i)^ for the Landau (black dots) and the 
covariant gauge (crosses) data. The residual gauge-dependence is of the order of 1.5% 



copies in the statistical sample of configurations 
generated and for the lattice used. 

To show our sensitiveness to the choice of the 
gauge, we compare the results obtained in the 
Landau and in the covariant gauge with the bare 
gauge parameter fixed to f = 8. 

The primary quantities measured on the lat- 
tice, Z^, A A and As, show a statistically signifi- 
cant gauge dependence (as you can see in Fig. 2 
of Ref. ^) while the renormalization constants 
Za, Zs shown in Fig. ||, obtained by computing 
the ratios as indicated in Eq. (Q), could have at 
most a residual gauge dependence (of the order of 
1.5%) which is not detectable within the statisti- 
cal precision of our simulation. For these quan- 
tities, the fluctuations of the simulation hide the 
weak gauge dependence that is expected in per- 
turbation theory from the next-to-leading order 
terms in a^. 

These results indicate an upper limit to the 
numerical troubles that can be expected in the 
RI/MOM renormalization constants because of 
the lattice gauge fixing. 
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